Equilibration of a spinless Luttinger liquid 
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We study how a Luttinger liquid of spinless particles in one dimension approaches thermal equi- 
librium. Full equilibration requires processes of backscattering of excitations which occur at energies 
of order of the bandwidth. Such processes are not accounted for by the Luttinger liquid theory. We 
treat the high-energy excitations as mobile impurities and derive an expression for the equilibration 
rate in terms of their spectrum. Our results apply at any interaction strength. 
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The concept of Luttinger liquid was proposed by 
Haldane as an effective low-energy description of one- 
dimensional systems of interacting fermions |l| or bosons 
[2| . The main feature of this theory is that regardless of 
the statistics of the particles, the low-energy excitations 
of the system are bosons. The latter propagate at a fixed 
velocity v in either left or right direction and have the 
meaning of the waves of particle density, analogous to 
phonons in solids. 

In its simplest form the Luttinger liquid is described 
by a Hamiltonian quadratic in boson variables, result- 
ing in excitations with infinite life time. Once excited, 
such a system will never reach thermal equilibrium. Ab- 
sence of equilibration is the physical reason [3j for per- 
fect quantization of conductance of a quantum wire con- 
nected to ideal leads, when the electronic system in the 
wire is treated as a Luttinger liquid [4|. Of course, real 
systems do equilibrate, possibly explaining the experi- 
mentally observed corrections to quantized conductance 
[5( . Equilibration of one-dimensional boson systems was 
recently studied in atomic traps [6(. 

A finite life time of excitations in the Luttinger liq- 
uid can be understood if small anharmonic corrections 
are added to the Hamiltonian. Such perturbations are 
irrelevant in the sense that their effect rapidly decreases 
as the temperature approaches zero. However, they are 
responsible for the interaction of bosonic excitations and 
therefore for their equilibration. Scattering of the excita- 
tions caused by the anharmonic coupling terms preserves 
not only their total energy but also momentum. Thus 
the resulting equilibrium distribution of the bosonic ex- 
citations 
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is controlled by two parameters, temperature T and ve- 
locity u. Here q is the wave vector of the excitation. 

It is important to note that translation invariance of 
the problem ensures conservation of the total momen- 
tum of the system, rather than that of its elementary 
excitations. This subtle distinction can be understood 



by considering the expression 

_ irhN T ^ . 
P=—-J+J2hq, 
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for the momentum of a Luttinger liquid [l|, [2J . Here b q 
is the boson annihilation operator, N is the total num- 
ber of particles, L is the system size. Periodic boundary 
conditions require that J be an even number if the un- 
derlying physical particles are bosons, while for fermions 
J + N must be even. The first term in Eq. @ accounts 
for the momentum associated with the motion of the sys- 
tem as a whole, which is possible even in the absence of 
excitations. 

Unless additional conservation laws are present, one 
should expect the existence of scattering processes which 
transfer momentum between the excitations and the sys- 
tem as a whole. The minimum momentum transfer 
Ap = 2irHN/L corresponds to J changing by 2. Be- 
cause the typical momentum of an excitation hq ~ T/v 
is small at T — > 0, such processes involve a large number 
of excitations. They are not included in the standard 
Luttinger liquid theory. Although their rate is small, 
these processes are required for the full equilibration of 
the Luttinger liquid. Physically one expects them to lead 
to relaxation of the velocity u of the gas of excitations in 
Eq. {!]) towards an equilibrium value Vd, 
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We limit our consideration to Galillean invariant systems 
of particles, whose mass is denoted by m. In this case 
the system must be at rest in a reference frame moving 
with the center of mass, and Vd = P/mN. The study of 
the relaxation time r is the main goal of this paper. 

We start by reviewing the simplest case of a system 
with Luttinger liquid behavior at low energies, namely, 
the weakly interacting Fermi gas. In fermionic Lut- 
tinger liquids the integer J can be interpreted [l|] as 
the difference of the numbers of right- and left-moving 
fermions, J = N R — N L . Clearly, the scattering pro- 
cess changing J by 2 involves backscattering of a fermion, 
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FIG. 1: (a) The simplest backscattering process in a weakly 
interacting Fermi gas involves three particles, including one 
near the bottom of the band, (b) Backscattering of a fermion 
at the bottom of the band can be interpreted as a hole exci- 
tation overcoming a barrier at Q — kp. 



AN R — —AN L = ±1. Because of conservation of energy 
and momentum, two-particle scattering in one dimen- 
sion results only in particles exchanging their momenta, 
and the distribution function remains unchanged. Thus 
the simplest scattering process involves three particles, 
see Fig. QJa) . Simultaneous conservation of momentum 
and energy requires involvement of hole states below the 
Fermi level. At T — > the most efficient process involves 
a hole near the bottom of the band, whose scattering 
is accompanied by creation and collapse of particle-hole 
pairs with energies of order T near the two Fermi points 
[7| . Since the occupation probability of a hole state near 
k = is exponentially small, one finds a small equilibra- 
tion rate t _1 oc e~ EF ^ T ', where Ep = h 2 k F /2m is the 
Fermi energy [8(. 

It is instructive to rephrase the above argument in the 
language of a hole excitation with wave vector Q and 
energy e(Q) — HvfQ(1 — Q/2kp) constructed by moving 
a fermion from state kp — Q to the Fermi level state kp. 
(Here vp is the Fermi velocity.) The hole is scattered 
off of particles near the Fermi level, with its momentum 
changing in steps of AQ ~ T/vp. Backscattering occurs 
when such a hole crosses the point Q = kp, Fig. [TJb). 

This picture can now be generalized to the case of ar- 
bitrary interaction strength. The particle hole pairs with 
momenta hq ~ T/v near the two Fermi points transform 
into the bosonic excitations in the Luttinger liquid [l|. 
On the other hand, the hole with the large wave vector 
Q ~ kp is not accounted for by the Luttinger liquid the- 
ory and should be treated as a mobile impurity [9|, |l(| . In 
the presence of interactions its energy e(Q) is defined as 
that of the lowest energy state of momentum hQ, mea- 
sured from the ground state. Throughout this paper 
we assume that e{Q) remains convex. Then the equi- 
libration rate shows activated temperature dependence 
t _1 oc e~ e ( kp >' T , where kp — ttuq is determined by the 
average particle density uq — N/L. 

To obtain a full expression for the equilibration rate, 
the distribution function of the holes should be consid- 
ered carefully. To first approximation it can be obtained 
by noticing that the holes are scattered by the bosonic 



excitations, distributed according to Eq. ([]}. These scat- 
tering events involve exchange of both energy and mo- 
mentum between the hole and the bosons, leading to the 
equilibrium distribution 



f(Q) 



^Q)/ T , Q < k F , 
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where e u (Q) = e(Q) — JiuQ. The apparent asymmetry 
between the cases of right- and left-moving holes, Q < kp 
and Q > kp, is caused by our convention to measure the 
momentum Q of the hole from the right Fermi point, 
k = +kp. 

The discontinuity of the hole distribution function (|4]) 
at Q = kp originates from the implicit assumption that 
the right- and left-moving holes are distinct particles. In 
reality, the backscattering processes shown in Fig. [1] con- 
vert right-moving holes into left-moving ones, thereby 
smearing the discontinuity of the distribution function 
f(Q). Because the hole moves in momentum space via 
random small steps of AQ ~ T/hv, this motion is dif- 
fusive. Such diffusion was considered previously for the 
cases of weakly-interacting [8| and strongly-interacting 
11] electrons. It is described by the Fokker-Planck equa- 
tion 



dtf = -d Q J, J = - 



B(Q) 



M 



T 



dc 



f, (5) 



where the expression for the probability current J as- 
sumes that the system as a whole is at rest, Vd = 0, 
and prime denotes the derivative with respect to Q. The 
diffusion constant in momentum space 



B(Q)=Y,[5Q] 2 W Q , Q+SQ 

8Q 



(6) 



is defined in terms of the rate Wq^q+sq of scattering 
events changing the wave vector of the hole from Q to 
Q + SQ. 

We now find a stationary solution of the Fokker-Planck 
equation with the boundary conditions ((31), which gives 
a uniform in Q-space probability current 



J = uB(kp) 



hkp 
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Here to obtain the expression for J to leading order in u 
we neglected the difference between e(Q) and e u (Q). 

A non-zero probability current J means that the holes 
backscatter at a rate JL/2tt. With each backscattering 
event transferring momentum Ap = 2Kkp from excita- 
tions to the motion of the system as a whole, we find 
P CK = —JLhkp/n. Comparing this result with the ex- 
pression P cx = (irLT 2 /3hv 3 )u for the total momentum 
of the excitations obtained using the distribution ([1]) , we 
find the relaxation law u = —u/t with the rate 
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Here A = e(pp), the effective mass of the hole to* = 
—H 2 /e"(kF), and the diffusion constant B = B(kp) re- 
mains to be determined. 

Following Refs. [9|, |l0|, we treat the hole in a Luttinger 
liquid as a mobile impurity. The Fokker-Planck equation 
for such an impurity was discussed in Ref. 12]. The 
parameter B was found to scale as 



B = X T b 



(9) 



at T — > 0. The approach of Ref. [12j does not allow 
for the determination of the coefficient %. The latter is 
controlled by the interactions between the physical par- 
ticles forming the Luttinger liquid. In the limit of strong 
Coulomb repulsion it was calculated in Ref. |ll|. A re- 
lated calculation was performed in the context of decay of 
dark solitons in weakly-interacting one-dimensional Bosc 
systems [131 ] - 

Our next goal is to obtain an exact expression for the 
coefficient \ m Eq. © for arbitrary interactions between 
the particles forming the Luttinger liquid. Microscop- 
ically the case of arbitrary interaction strength can be 
approached only for integrable systems, where an infi- 
nite number of conservation laws allows one to diagonal- 
ize the Hamiltonian exactly. However, the same conser- 
vation laws ensure that the excitations have infinite life 
times and B — 0. We thus develop a phenomenological 
theory and express B in terms of hole spectrum e(Q). 

We describe the system in terms of the displacement 
u(y) of a small element of the liquid from its refer- 
ence position y in a state of uniform particle density 
no, and the conjugate momentum density p(y) such that 
[u(y),p(y')] — ihS(y — y'). In the absence of the hole 
excitations the Hamiltonian of the liquid can be written 
as 



Hr = 



ImriQ 



n U(n) 



dy, 



(10) 



where U(n) is the internal energy per particle, deter- 
mined by the fluctuating density n{y) = Uo/[l + u'{yj\. 
Expanding (jlQI) up to the third order in small deforma- 
tion vl one finds 



-f(X 



hl= j u^ + ^ u ~ au ) y - (n) 

Here the sound velocity v = [(2n U' + nlU'^/m] 1 / 2 and 
a = n\XJ' + n%U" + n^U'"/6. The quadratic part of 
Eq. (TTTJ) is the Hamiltonian of the Luttinger liquid, which 
can be brought to the form ^2hv\q\b^b q by introducing 
the boson operators b q via the standard procedure 



«(v) = E 



2mnoLv\q 



(6,e*w + 6je- 4w ), (12a) 



^) = -E^^(v^-^-). ( i2b) 



The presence of a hole excitation at the point in the liq- 
uid with reference position Y is accounted for by adding 
a term Hh = e(Q) = e(— idy) to the Hamiltonian (|11|) . 
Since our goal is to evaluate B = B(fcp-), we assume that 
Q is near kp — ttHq and use the expansion 



H h = A(nQO) - 7^ — (-Wr - ™o) 

2m* 



(13) 



It is worth mentioning that our Hamiltonian is written in 
terms of the Lagrangian variable u(y) defined as function 
of reference position y, rather than Eulerian variable n(x) 
at the physical position x — y+u(y). The two approaches 
are, of course, equivalent and lead to the same results 
14] . Although the use of Eulerian variables is more com- 
mon in the Luttinger liquid theory, our method has the 
advantage of more simply accounting for the Galilean 
invariance of the problem. In addition, since Y is the 
position of the hole in the reference state of uniform den- 
sity no, the maximum of e(Q) is located at Q = ttuq, 
regardless of the physical density n. On the other hand, 
the maximum value A is a function of n = no/(l + u'). 
This dependence gives rise to interaction of the hole with 
the Luttinger liquid. Expanding (|13|) to second order in 
u' , we obtain 

H h = ^(3 1 u'(Y)+p 2 [u'(Y)} 2 -^- ; (-tdY-7Tn ) 2 , (14) 

where j3\ = tojA', fa — «oA' + nQA"/2, and we omitted 
the constant A(no). 

In order to find the diffusion constant in momentum 
space B(kp), Eq. ((6]), we evaluate the scattering rate 
Wq^q+sq- The momentum of the hole changes as it inter- 
acts with the bosonic excitations, see Eqs. ([T4|) and (fT2|) . 
The processes involving one boson cannot simultaneously 
conserve both energy and momentum of the system. The 
simplest allowed process for a hole near Q — ttuq involves 
absorption of a boson q\ and simultaneous emission of a 
boson Q2 such that q 2 ~ —qi [3, LL2, LUa] • The scattering 
rate is then found from the Fermi golden rule expression 

W Q,Q+SQ = -r- z2 1*91-92 1 2 N qi( N Q2 +^)Sq 1 -q 3 ,SQ 

x«5(e(Q) - e(Q + SQ) + hv\ qi \ - hv\q 2 \). 

The matrix element t qi _ q2 accounts for all processes that 
destroy boson q\ and create boson q 2 . For example, a 
contribution proportional to fabl b qi is found in the sec- 
ond term in Eq. dT4l) . Identical scattering processes can 
be obtained in the second-order perturbation theory with 
amplitudes proportional to 01 or afti . The calculation is 
simplified considerably by applying to the Hamiltonian 
the unitary transformation U'(Hl + Hh)U with 



U = exp 



tPi 



hmnov 2 



P(y)dy 
J 



(15) 



This removes the —fau'(Y) term in (fT4)) and generates a 
correction to fa proportional to a fa . In addition, a new 
term p 2 (Y) is generated with the coefficient proportional 
to fay. Both the [u'iY)] 2 and p 2 (Y) terms contain con- 
tributions of the form 6j b qi and give rise to the matrix 
element 



"9i,92 



TV I gi 92 1 /„ 3a/3i 
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ft 
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(16) 



As a result, we recover the temperature dependence Q 
with the coefficient x given by 



4ir 3 n 2 



a// 2w'*/ A ' 2 

X = s — - — I A" A' H ^ 



(17) 



where prime denotes the derivative with respect to the 
particle density uq. 

The above result completes our evaluation of the re- 
laxation rate of a Luttinger liquid, given by Eqs. ([8]). 
(O, and (fTT|) . The rate has activated temperature de- 
pendence with both the activation temperature A and 
the prefactor determined by the spectrum of holes e(Q). 
Although our result is applicable at any interaction 
strength, the spectrum e(Q) is known analytically only 
in a few special cases. 

For non-interacting spinless fermions A is given by 
the Fermi energy (irhno) 2 /2m, v is the Fermi velocity 
Trhno/m, and m* = m. This results in x = 0, as there is 
no scattering of holes in the absence of interactions. In 
the limit of weak interactions, the spectrum e(Q) should 
be evaluated up to second order in interaction strength. 
This gives rise to a result [14J for \ consistent with the 
rather complicated expression for the three-particle scat- 
tering amplitude [7] that controls the scattering of holes, 
Fig. [lja). In the limit of strong long-range repulsion, the 
system forms a Wigner crystal, and the hole spectrum 
coincides with that of phonons in the crystal. We have 
verified that in this regime our expression (|17|) recovers 
the results of Ref. [llj. We have also found that in the 
case of weakly interacting bosons Eq. (fl7|) is consistent 
with the expression for the mobility of the so-called dark 
soliton 11311 . 



At arbitrary interaction strength the spectrum of holes 
is known only for integrable models. As we already men- 
tioned, integrability means absence of scattering of ex- 
citations, B = 0. We have verified that our expression 
(fTTf vanishes for the Calogero-Sutherland model of par- 
ticles with inverse-square repulsion [15[ , and for the Lieb- 
Liniger model of bosons with point-like repulsion [16j . 
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